Two-dimensional Dirac fermions with random axial-vector potential 
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A Dirac fermion model with random axial-vector potential is proposed. At a special strength 
of randomness, the symmetry of the action is enhanced, which is due to the gauge symmetry a 
la Nishimori. Some exact scaling exponents of single-particle Green functions are computed. The 
relationship with the XY gauge glass model is discussed. 
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Random Dirac fermions have attracted much current 
interest in condensed matter physics. They have actually 
intimate relationship with integer quantum Hall (IQH) 
transitions dirty d-wave superconductors l^,!^!, etc. 

Critical theory of IQH transitions is believed to be a 
strong coupling fixed point of the Dirac fermion with 
random mass, vector and scalar potentials though the 
fixed point is still missing. Near the zero energy quasi- 
particles of disordered d-wave superconductors can be de- 
scribed by the Dirac fermion, giving four kinds of novel 
universality classes 1^,0. Including the conventional or- 
thogonal, unitary and symplectic classes and correspond- 
ing chiral classes, ten universality classes have been es- 
tablished for disordered systems 

On the other hand, it has been reported rec^cntly that 

~ ~0). One 
llfof the 



several other symmetry classes should exist 
example is realized |^ on the Nishimori-linc 
random-bond Ising model [[l2[. The random-bond Ising 
model has 0(2n) symmetry ||13| , p^ , belonging to the uni- 
versality class D (the class of superconductors with 
broken spin rotation and time reversal symmetry) in 
Zirnbauer's classification Q. Although the enhancement 
of the symmetry on the Nishimori-linc has already been 
discussed in Gruzberg et. al. have exphcitly shown 
that the symmetry is enhanced to 0(2n -I- 1) due to the 
gauge symmetry This class is not included in Zirn- 
bauer's classification. 

Among these developments, the Dirac fermion with 
random vector potential only is one of simplest but non- 
trivial models, still providing hot topics. It gives mul- 
tifractal scaling exponents of local composite operators 
, exact zero-energy wave function for any realization 
of disorder |p],p|,p^ , rephca symmetry breaking |]l7| , |l8| , 
etc. 

In this paper, we study a Dirac fermion in two dimen- 
sions including random axial-vector potential. We show 
that the symmetry of the model is enhanced at a strong 
disorder strength, which is analogous to the Nishimori- 
line of statistical models. Similarly to the conventional 
model with random vector potential, this model can be 
solved exactly, giving us some exact scaling dimensions. 
It turns out that this model corresponds to the spin wave 
model for the XY gauge glass model. 

The system we will study is 



Ha = -tphAip, liA = 7m {-idf, - y^jsAfj 



(1) 



where 71, 72, and 75 denote the three Pauli matrices, 
and Af^ is a random axial-vector potential. Being de- 
fined in 2D Euclidean space, this Hamiltonian is non- 
hermitian due to 75 coupled with the vector potential. 
Therefore, the model (H 
Hatano-Nelson model [1 



may have a relationship with 
. However, ensemble- average 
yields the square of the antihermitian term, and the La- 
grangian becomes hermitian. It should also be mentioned 
that without mass, the Hamiltonian has chiral symmetry 
75 + ^A75 = 0, telling that the model belongs to the 
class ^III in Zirnbauer's classification. 

In the previous work the same model but with 

usual vector potential, hy — '^f_i. ~ \/9^fj) '^^^ 

studied. In the present case, the coupling with 75 may 
imply that the right- and left-moving fermions have op- 
posite charges. In what follows, the present model (|^) is 
sometimes referred to as A-model, whereas the conven- 
tional one as V-model. The probability distribution of 
the axial-vector potential is assumed to be of Gaussian 
type 



P\A] = —e-^I''' 
^ ^ Na 



(2) 



where the normalization factor is given by Na — 

At a given Matsubara frequency lu, 
the partition function of the model with quenched disor- 
der is given by 



Z 



(3) 



The Lagrangian is invariant under the following chiral 
gauge transformation. 



A^ ^ A^ + df,e{x). 



(4) 



It is well-known that this transformation gives rise to a 
nontrivial Jacobian, but as we shall see later, it is irrele- 
vant to our problem. The w-term serves as a symmetry- 
breaking term. 

Although it is possible to take quenched average of 
the partition function without resort to the replica 
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trick, we will first apply it in order to explain the 
enhancement of the symmetry at a special coupling 
constant. Integration over disorder for replicated 
partition function yields an effective Lagrangian, 

Ca = V'" {-Jtiid^ + iuj) I/'" - f (V'"757mV'°) , where a = 
l,...,7i. Due to the identity ^^.js — —i^^ivlv, we 

have — I {4>°'lblii'^°') — +f (V'^Tm^") • Compared 
with the effective Lagrangian of the V-model, Cy = 

{—Ifiid/j. + iio) V''^ — f (V'°7^V''') J the present A-model 
turns out to be equivalent to the V-model but with nega- 
tive coupling constant g. Namely, the partition function 
Z^{g) of the A-model with the coupling constant g > 
is equivalent to that of the V-model Zy{—g). Therefore, 
some properties, e.g., the scaling dimensions of some lo- 
cal operators for the A-model can be obtained directly 
from the V-model by the use of the formal replacement 
g —g, as we shall see momentarily. 

However, the A-model has some peculiar properties. 
One of them is the enhancement of the symmetry at 
a special coupling constant g — t:. The Lagrangian in 
Eq. (H) has global \J{ri)x \J{ri) symmetry for any g 
when uj = 0. Explicitly, it is readily seen that for any 
realization of the axial-vector potential it is invariant 
under the transformation {UiV+ + U2'P-)ip and 

^p ■ip{Ulr-+UlV+), where V± = (l±75)/2 and Ui is 
a n X n unitary matrix which acts on the replica space. 
Since the Lagrangian is invariant under local chiral gauge 
transformation, we have 



to n, whereas a from 1 to n. When g = n (and 
Lu = 0), the ensemble-averaged Lagrangian can be writ- 
ten in a symmetric way together with the 0th fermion 

as £ = -V>7m«'9/xV'" + i (V'"7MV'")^ telling that the 
symmetry is enhanced to U(n -I- l)x U(rt -I- 1). This en- 
hancement of the symmetry for the Dirac fermion with 
random axial-vector potential is quite analogous to that 
for the Majorana or Dirac fermion with random mass 
(the random-bond Ising model). Interestingly, such an 
enhancement does not occur in the case with conven- 
tional vector potential. 

What happens at this special coupling constant? For 
the conventional model hy, the exact scaling properties 
have been obtained by various methods Since the 

present Hamiltonian Ha is nonhcrmitian, we exactly solve 
the model by taking the quenched average directly, with- 
out using the replica trick or the SUSY technique. One 
of most interesting properties in disordered systems is 
the multifractality of the scaling exponents. The Dirac 
fermion with random vector potential is a typical exam- 
ple having such exponents. By the use of the single- 
particle Green function. 



Gtj{x,y,iu}) = {'^/J^{x)^pJy)), 



where 



-J d^xC 



(8) 



(9) 



Z"= / V^pVipVAe 



-\ J d^^Al - J <fxC 
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(5) 



where C — {id^ + -^75^^) ip"', V is the "vol- 

ume" of the gauge space V — J "DO, and r(^is a Ja- 
cobian due to chiral gauge transformation ||2^. This 
equation is based on the same argument as Nishimori 
applied to the random-bond Ising model Since the 
Jacobian is multiplied by a factor n, it should vanish in 
the replica limit n — > 0. This is also justified below in 
the direct computation of the exact scaling dimensions 
of the single-particle Green functions. Now using the 
fermionization of the gauge variable via fermion-boson 
correspondence in 2D, 

i(a^^)2 ^-V^V^mV'", ^e^.a,0 ^ V^V^o, (6) 



we have 



where the Lagrangian now includes an additional 
0th fermion £ — —{p^j^idfj^tlj" — V5(V''^7m75V''' + 
./ ^tp'~''yfj,'y5ip^)Af^. The Greek superscript a runs from 



ensemble-averaged density of state (DOS) 
and a field-theoretical analogue of the in- 
verse participation ratios (IPR) P^'^'^{E) = 
En |*„(a:)P'=i£; - i?„)/E„ \^nix)\^SiE ~ i?„) can be 
computed as |gl| 



piE) 



— lim Im trG(x, x^iuj — E), 



P^^\E)p{E) = -i- \iuiLo^'^[luiirG[x,x,E ~iuj)Y, (10) 

where ^n(a:) is a two component spinor wave function 
of Tith eigenstate, tr is a trace for spinor indices, and 
Cfc = iT{2k - 3)!!/(2n - 2)!!. It is useful to introduce 
the chiral basis, ip± — V±ip and i/'l — V'T^ip, which was 
actually used by Mudry et. al. ||]. In what follows, 
we will take similar notations to this reference. Since 
■0^/; = -|- tpl^^^, the kth power of the Green func- 

tion can be written as 



[G{x,x, E - cx ( n [^l"C(2:) + ijl^r+i^)] ). (11) 

a=l 

Here we have introduced species a to calculate the A:th 
power of the Green function, although we will not use the 
replica trick. Expanding the product, we have various op- 
erators with various scaling dimensions. Most dominant 
operators depends on whether k is even or odd. Define 
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[fe/2] 

0(fe)(x) = II Vi^'-V-'-Vi'V^lx), (12) 

a=l 

then an example of most dominant operator is Oe'^'' = 
OC^) for even k and oi'''^ = 0('=)V+ V'- for odd k. 
In the chiral basis the Lagrangian is converted into 

C = ^ljl^{2td, - ^A+Wl + yjl^2td, + ^A_)r-, (13) 

where z = x + iy, z ~ x — iy, A± ~ Ax ± iAy^ 
and a — l,2,---,k. The different sign for A± tells 
that the right- and left-movers have opposite charges, 
as mentioned previously. In 2D the gauge fields can 
be decomposed into two independent components — 
^ni^d^rj + d^S,, or in the chiral basis, A^ = —2idz{r/ + iQ, 
and A^ — 2idz{ri — The probability distribution (g) 
is written now in terms of these fields as 

where ^i,r) is a normalization factor iV^ = 
/P^e"*/'^'"'^^''^^' and similar for Nr,. Via the 
gauge transformation ?/;^ gVaiv±iS.) ^i^^ and ip^t ^ 
the axial-vector potential completely dis- 
appears in the Lagrangian, 

V^"2i9ji/;^-K7/'!_''2ia^V-- (15) 

The chiral gauge transformationyields a nontrivial Jaco- 
bian VipVip V^jV-ipe'-''^'^^'^ [|o|. However, this factor 
appears in the denominator Z as well as in the numer- 
ator in Eq. (^), and hence cancels out. Finally, Oe'^^ is 
gauge-invariant while Oo*'^ is not. Actually, they obey 
the transformation laws 

oik) _^ Q(k)^ Q{k) _^ e-2iy5?o(fc)^ (16) 
By the use of the free action ( |l5| ) after the gauge 

ik'] 

transformation, the correlation function of Oe,o can 
be evaluated as follows. The two point correlators of 
the free fermi fields are {ip'\,{z)il}]^ (0)) ~ (S"" z~^ and 
(■0^(z)'0l°'(O)) S^^'z-^. Therefore, we have 

mx)oim = \x\-^\ 

{Oo{x)Ol{0)) =e-2»v^«(-)e2»Vff«(o)|2;|-2fc^ ^^7) 

Note that these do not depend on -q and therefore the 
change of the scaling dimension due to disorder is in- 
volved with ^ only. The average-over ^ is easily taken 

ik') 

and the scaling dimension of Oe,o finally reads Afe — k 
for even k whereas Ai. — k+^ for odd k, and this is itself 
nothing but the dominant scaling dimension of the kth 
power of the single-particle Green function (trG)'^. Some 
comments may be in order. In the conventional V-model 
the random field rj plays a role in the change of the di- 
mensions due to disorder. In the present case, the role 



of r] and ^ is exchanged by 75 in Eq. (P, and ^ causes 
the change of the scaling dimensions. Furthermore, the 
field 77 gives in general negative scaling dimensions for 
the V-model, whereas ^ gives positive dimensions for the 
present model. 

Using the dominant scaling dimensions of the Green 
functions obtained so far, we can calculate the scaling 
properties of the DOS and IPR. Since uj couples with '0V' 
in the Lagrangian, the dimension of uj reads z = 1 — ^ 
||l|,||. Therefore, we expect p{uj) ~ a;^^"^^/^. This is just 
the same formula for the V-model but with negative g. It 
is also readily seen that the random axial- vector potential 
is a marginal perturbation and the theory moves along 
a critical line. However, starting from the free fermion 
point, the A-model moves, as g increses, to the opposite 
direction to which the V-model moves. Therefore, the 
line on which the A-model lies is not reached by the con- 
ventional random vector potential model. In this sense, 
the present model is complementary to the full critical 
line of U(l)xU(l) symmetry. The dimension of p is an 
increasing function of g, so that the DOS is suppressed 
around the zero energy if the random axial-vector poten- 
tial is switched on. At the special point g = tt where 
the symmetry is enhanced, the exponent becomes infin- 
ity. Therefore, it is likely that the theory is well-defined 
only for g < n, and the enhancement of the symmetry is 
a signal that the theory reached at a singularity. 

The scaling behavior of the IPR is obtained in a sim- 
ilar way: P'*^' ^ cj'^'--/^, where — —2 for even k and 
= for odd k. As stressed by Mudry et. al. the 
IPR defined by Eq. (|l^) is not necessarily coincides with 
that of the original definition pf*^^ = |^'„(x)/A/'p'^, where 
the normalization A/" is defined by TV^ = J (Px\'i/o{x)\'^ . 
Actually, the zero energy wave function can be obtained 
exactly as 

^oix) = e^^^^-^^t^^^Vo, (18) 

for the present A-model, where (po is a constant spinor, 
while ^'0(2;) = e''^^)^'^~*^(^Vo for the conventional V- 
model. The IPR of the latter wave function has been 
studied by various methods, e.g., repHca method [Q, 
SUSY method §, the method using the equivalence to 
the random energy model ]l^, etc. They have all given 
the same scaling exponent for small g. Eq. ( [T^ ) for the 
present model also gives |^'o(a;)p oc e^'''^-'. Since the 
present A-model has the same probability distribution 
( p^ ) as the conventional V-model has, we expect that 
the the A-model should have the same scaling dimension 
Tk — {2 — ^k){k — 1) as the V-model has, which is quite 
different from the conjectured based on ([ic|). 

So far we have studied some peculiar properties of the 
Dirac fermion with random axial- vector potential. As al- 
ready mentioned, the Hamiltonian is nonhermitian, be- 
cause it is defined in the Euclidean space. Then, is this 
model unrealistic? To address the question, we next ex- 
amine the bosonized form of the A-model, as has been 
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done for the V-model by Bernard g| . By the use of the 
correspondence , the Hamiltonian ([^) is converted into 



n 



(19) 



This is a low-temperature effective Hamiltonian of the 
XY gauge glass model ||ll|,^,^ whose Hamiltonian is 
defined on a 2D square lattice, 



[3H cos((/), - + Xij). 



(20) 



where Xij is a random gauge variable with a probability 
distribution 



(21) 



Provided that K and Kp are large, we actually reach 
the Hamiltonian (^9|) with g = ^ and the probability 

distribution (|^) as well by expanding the cosine terms 
up to second order, defining (f>i — (j)- 



a{)d^(j){x) and 

Xij ~ o,oA^{x) with the lattice constant oq, and rescal- 
ing the fields (j) <j>/VK and A^J ^JlTp. This 

is the so-called spin-wave approximation. The spin-wave 
Hamiltonian ( p^ is expected to be unstable against vor- 
tex excitations if K or Kp becomes small. 

On the other hand, it has been shown that the model 
( pO| ) with (|2l| ) has a Nishimori-line Kp = K, on which 
some exact results can be obtained in a similar way as 
the random-bond Ising model ||ll|,^. The equivalence 
of the condition Kp = K to g = tt implies that the 
symmetry enhancement of the Dirac fermion with ran- 
dom axial-vector potential reflects the gauge symmetry 
on the Nishimori-line of the gauge glass model. As to the 
symmetry of the model, it is difficult to read its enhance- 
ment from the spin- wave Hamiltonian (19). Even if we 
apply the replica method, the symmetry always remains 
U(l). This reminds us of the case of the random-bond 
Ising model: We only come accross the continuous 0(n) 
symmetry when the model has been described by the Ma- 
jorana fermion via the Jordan- Wigner transformation. In 
the present case, the enhanced symmetry also becomes 
manifest only in the fermion description. However, the 
Lagrangian discussed so far is just for the spin- wave part 
of the gauge glass model and it is quite necessary to in- 
clude vortices to fully describe the lattice model (20). 
Then, such theories should be well-defined field theories 
with manifest enhanced symmetry on the Nishimori-line. 
It is an interesting issue to derive and study such field 
theories. 

Finally, let us discuss the symmetry for the dual the- 
ory of the gauge glass model (pO|) |Q . The periodicity of 
the Hamiltonian can be described by '^"^'^'^^^'^^^^^^^ ^ 
Em=-oo e'^e--T('^--'^^+X'^-2^™). By ^j^g of ^j^g pQjg_ 
son summation formula 124], we reach the following 
Hamiltonian, Hd — ^{dfjOj^ + ief^^Af^d^O — 2?/cos27r6', 
where is a dual field and y is a fugacity. This is the sine- 
Gordon model coupled with random vector potential A^. 



Therefore, it turns out that there is a duality between 
the random axial-vector potential and conventional vec- 
tor potential. The former model yields a symmetry en- 
hancement whereas the latter does not. Recently, the 
replica symmetry-breaking of the V-model has been sug- 
gested ||l3,|l^ , which might be understood via this duality 
relation. 

In summary, we have studied a Dirac fermion with ran- 
dom axial-vector potential. This model has an enhanced 
symmetry at a special strength of randomness, which re- 
minds us of the Nishimori-line of the statistical models. 
Indeed, the model is equivalent to the spin wave Hamilto- 
nian of the XY gauge glass model. It turns out that this 
model moves along a critical line for increasing g up to 
TT, and some exact scaling exponents have been obtained. 
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